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Abstract. We study the global dynamics of the universe within the framework of the Interacting Dark Matter (IDM) scenario. 
Assuming that the dark matter obeys the collisional Boltzmann equation, we can derive analytical solutions of the global density 
evolution, which can accommodate an accelerated expansion, equivalent to either the quintessence or the standard A models, 
with the present time located after the inflection point. This is possible if there is a disequilibrium between the DM particle 
creation and annihilation processes with the former process dominating, which creates an effective source term with negative 
pressure. Comparing the predicted Hubble expansion of one of the IDM models (the simplest) with observational data we find 
that the effective annihilation term is quite small, as suggested by a variety of other recent experiments. 
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1. Introduction 

Over the past decade the analysis of high quality cosmolog- 
ical data (supernovae type la, CMB, galaxy clustering, etc.) 
have suggested that we live in a fiat and accelerating universe, 
which contains cold dark matter to explain clustering and an 
extra component with negative pressure, the vacuum energy 
(or in a more general setting the dark energy), to explain the 
observed accelerated cosmic expansion (Spergel et al. 2007; 
Davis et al. 2007; Kowalski et al. 2008; Komatsu et al. 2009 
and references therein). Due to the absence of a physically 
well-motivated fundamental theory, there have been many the- 
oretical speculations regarding the nature of the above exotic 
dark energy (DE) among which a cosmological constant, scalar 
or vector fields (see Weinberg 1989; Wetterich 1995; Caldwell, 
Dave & Steinhardt 1998; Brax & Martin 1999; Peebles & Ratra 
2003; Perivolaropoulos 2003; Brookfield et al. 2006; Boehmer 
& Harko 2007 and references therein). 

Most of the recent papers in this kind of studies are based 
on the assumption that the DE evolves independently of the 
dark matter (DM). The unknown nature of both DM and DE 
implies that we can not preclude future surprises regarding 
the interactions in the dark sector. This is very important be- 
cause interactions between the DM and quintessence could pro- 
vide possible solutions to the cosmological coincidence prob- 
lem (Grande, Pelinson & Sola 2009). Recently, several pa- 
pers have been published in this area (eg., Amendola et al. 
2003; Cai & Wang 2005; Binder & Kremer 2006; Campo et al. 
2006; Wang, Lin & Abdalla 2006; Das, Corasaniti, & Khoury 
2006; Olivares, Atrio-Barandela & Pavon; He & Wang 2008, 



and references therein) proposing that the DE and DM could 
be coupled, assuming also that there is only one type of non- 
interacting DM. 

However, there are other possibilities: (a) It is plausible that 
the dark matter is self-interacting (IDM) [Spergel & Steinhardt 
2000], a possibility that has been proposed to solve discrep- 
ancies between theoretical predictions and astrophysical obser- 
vations, among which the less cuspy halo profiles, predicted 
by the IDM model, allowing for the observed gamma-ray and 
microwave emission from the center of our galaxy (Flores 
& Primack 1994; Moore et al. 1999; Hooper, Finkbeiner & 
Dobler 2007; Regis & Ullio 2008 and references therein) and 
the discrepancy between the predicted optical depth, r, from 
the Gunn-Peterson test in the spectra of high-z QSOs and the 
WMAP-based value (eg., Mapelli, Ferrara & Pierpaoli 2006; 
Belikov & Hooper 2009; Cirelli, Iocco & Panci 2009 and ref- 
erences therein). It has also been shown that some dark matter 
interactions could provide an accelerated expansion phase of 
the Universe (Zimdahl et al. 2001; Balakin et al. 2003; Lima, 
Silva & Santos 2008), (b) The DM could potentially contain 
more than one particle species, for example a mixture of cold 
and warm or hot dark matter (Farrar & Peebles 2004; Gubser & 
Peebles 2004), with or without inter-component interactions. 

In this work we are not concerned with the viability of the 
different such possibilities, nor with the properties of interact- 
ing DM models. The single aim of this work is to investigate 
whether there are repercussions of DM self-interactions for the 
global dynamics of the universe and specifically whether such 
models can yield an accelerated phase of the cosmic expansion, 
without the need of dark energy. Note that we do not "design" 
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the fluid interactions to produce the desired accelerated cosmic 
evolution, as in some previous works (eg., Balakin et al. 2003), 
but rather we investigate under which circumstances the ana- 
lytical solution space of the collisional Boltzmann equation, in 
the expanding Universe, allows for a late accelerated phase of 
the Universe. 



the collisional Boltzmann equation reduces to the usual fluid 
equation: p + 3H(p + P) = 0, Inserting eqs.([3j and (0 into 
eq.(f2| we now obtain: 
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2. Collisional Boltzmann Equation in the 
Expanding Universe 

It is well established that the global dynamics of a homoge- 
neous, isotropic and flat Universe is given by the Friedmann 
equation: 



(£)-=?■ 
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with p the total energy-density of the cosmic fluid, containing 
(in the matter dominated epoch) dark matter, baryons and any 
type of exotic energy. Differentiating the above, we derive the 
second Friedmann equation, given by: 



a AnG ( p\ 
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As we have mentioned in the introduction, the dark mat- 
ter is usually considered to contain only one type of particle 
that is stable and neutral. In this work we investigate, using the 
Boltzmann formulation, the cosmological potential of a sce- 
nario in which the dominant "cosmic" fluid does not contain 
dark energy, is not perfect and at the same time it is not in 
equilibrium^. Although our approach is phenomenological, we 
will briefly review a variety of physically motivated dark matter 
self-interaction models, which have appeared in the literature. 

The time evolution of the total density of the cosmic fluid 
is described by the collisional Boltzmann equation: 



■£ + 3H(t)p + k P 2 -*¥ = 0, 



(3) 



where H(t) = a/a is the Hubble function, *F is the rate of cre- 
ation of the DM particle pairs and k(> 0) is given by: 



(cru) 
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(4) 



where cr is the cross-section for annihilation, u is the mean par- 
ticle velocity, and M x is the mass of the DM particle. 

Note that, in the context of a spatially flat FLRW cosmol- 
ogy, for an effective pressure term of: 
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1 Initially, the total energy density is p = piom + Pr- We consider 
that the self interacting dark matter does not interact significantly with 
the background radiation, and thus in the matter dominated epoch, ra- 
diation is irrelevant to the global dynamics (due to the well-known de- 
pendence: p r oc a~ 4 ). Therefore, taking the above considerations into 
account and assuming that there are no residual radiation products of 
the DM interactions (otherwise see Appendix A), we conclude that in 
the matter dominated era the total cosmic dark-matter density reduces 
to that of the IDM density (p s Pidm) which obeys the collisional 
Boltzmann equation (see eq.O. 



Obviously, a negative pressure (whichever its cause) can ef- 
fectively act as a repulsive force possibly providing a cosmic 
acceleration. 

In this paper we investigate the effects of DM self- 
interactions to the global dynamics of the Universe and under 
which circumstances they can produce a negative pressure and 
thus provide an alternative to the usual dark energy. It is well 
known that negative pressure implies tension rather than com- 
pression, an impossibility for ideal gases, but not so for some 
physical systems which depart from thermodynamic equilib- 
rium (Landau & Lifshitz 1985). 

The particle annihilation regime has been described by 
Weinberg (2008), using the collisional Boltzmann formula- 
tion, in which the physical properties of the DM interactions 
are related to massive particles (still being present) which if 
they carry a conserved additive or multiplicative quantum num- 
ber, would imply that some particles must be left over after 
all the antiparticles have annihilated (Weinberg calls them L- 
particles). The L-particles may annihilate to other particles, 
which during the period of annihilation can be assumed to be in 
thermal and chemical equilibrium (see Weinberg 2008). Such a 
DM self-interacting model has repercussions to the global dy- 
namics of the Universe (see our Case 2 below). 

The corresponding effects to the global dynamics of the 
particle creation regime, providing an effective negative pres- 
sure, has also been investigated by a number of authors (eg., 
Prigogine et al. 1989; Lima et al. 2008 and references therein). 

Generally, in the framework of a Boltzmann formalism, a 
negative pressure could indeed be the outcome of dark mat- 
ter self-interactions, as suggested in Zimdahl et al. (2001) 
and Balakin et al. (2003), if an "antifrictional" force is self- 
consistently exerted on the particles of the cosmic fluid. This 
possible alternative to dark energy has the caveat of its un- 
known exact nature, which however is also the case for all dark 
energy models. Other sources of negative pressure have also 
been proposed, among which gravitational matter "creation" 
processes (Zeldovich 1970), viewed through non-equilibrium 
thermodynamics (Prigogine et al. 1989) or even the C-field of 
Hoyle & Narlikar (1966). The effects of the former proposal 
(gravitational matter creation) on the global dynamics of the 
Universe have been investigated, under the assumption that the 
particles created are non-interacting (Lima et al. 2008). The 
merit of all these alternative models is that they unify the dark 
sector (dark energy and dark matter), since just a single dark 
component (the dark matter) needs to be introduced in the cos- 
mic fluid. 

In what follows we present, in a unified manner, the out- 
come for the global dynamics of the Universe of different type 
of dark matter self-interactions, using the Boltzmann formula- 
tion in the matter dominated era. 
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3. The Cosmic Density Evolution for different DM 
interactions 

We proceed to analytically solve eq. (f3]l. We change variables 
from t to a and thus eq.© can be written: 



^ = f(a)p 2 + g(a)p + R(a) 
da 



where 



/(«) 
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Within this framework, based on eqs.(|5j|7]i and ©, we can dis- 
tinguish four possible DM self-interacting cases. 

Case LP- 0: If the DM is collisionless or the collisional 
annihilation and pair creation processes are in equilibrium (ie., 
¥ = Kp 2 ), the corresponding solution of the above differential 
equation is p oc a~ 3 (where a is the scale factor of the universe), 
and thus we obtain, as we should, the dynamics of the Einstein 
de-Sitter model, with H(t) - 2/3t. 

Case 2: P - Kp 2 /3H: If we assume that in the matter era the 
particle creation term is negligible, *F = [R(a) - 0], (the case 
discussed in Weinberg 2008), then the corresponding pressure 
becomes positive. It is clear that eq. (0 becomes a Bernoulli 
equation, the general solution of which provides the evolution 
of the global energy-density, which is that corresponding to the 
IDM ansanz: 



p(a) 



C 2 - j" x- 3 f(x)dx C 2 + k£ a~ 3 (t)dt 



(9) 



Prior to the present epoch we have that p(a) oc a~ 3 , while at 
late enough times (a » 1) the above integral converges, which 
implies that the corresponding global density tends to evolve 
again as the usual dark matter (see Weinberg 2008), with 
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where to is the present age of the Universe. The latter analysis, 
relevant to the usual weakly interacting massive particle case - 
Weinberg (2008), leads to the conclusion that the annihilation 
term has no effect resembling that of dark energy, but does af- 
fect the evolution of the self interacting DM component, with 
the integral in the denominator rapidly converging to a constant 
(which does depend on the annihilation cross-section). 

Case 3: P - (Kp 2 - W)/3H: For the case of a non-perfect 
DM fluid (ie., having up to the present time a disequilibrium 
between the annihilation and particle pair creation processes) 
we can either have a positive or a negative effective pressure 
term. Although the latter situation may or may not appear plau- 
sible, even the remote such possibility, ie., the case for which 
the DM particle creation term is larger than the annihilation 
term (Kp 2 - *P < 0), is of particular interest for its repercus- 
sions on the global dynamics of the Universe (see for example 
Zimdahl et al. 2001; Balakin et al. 2003). 

It is interesting to note that this case can be viewed as 
a generalization of the gravitational matter creation model of 
Prigogine et al. (1989) [see also Lima et al. 2008 and references 



therein] in which annihilation processes are also included, al- 
though the matter creation component dominates over annihi- 
lations. In such a scenario, as well as in any interacting dark- 
matter model with a left-over residual radiation, a possible con- 
tribution from the radiation products to the global dynamics is 
negligible, as we show in appendix A. 

In general, for k + and *P + it is not an easy task to 
solve analytically eq. (|7]), which is a Riccati equation, due to 
the fact that it is a non-linear differential equation. However, 
eq.© could be fully solvable if (and only if) a particular so- 
lution is known. Indeed, we find that for some special cases 
regarding the functional form of the interactive term, such as 
*P = v P(a, H), we can derive analytical solutions. We have iden- 
tified two functional forms for which we can solve the previous 
differential equation analytically, only one of which is of inter- 
est since it provides a oc gt 3 dependence of the scale factor (see 
appendix B). This is: 



W(a) = aH(a)R(a) = C\(m + 3)a"'H(a) + K Cia 



2 „2m 



(ID 



Although, the above functional form was not motivated by 
some physical theory, but rather phenomenologically by the 
fact that it provides analytical solutions to the Boltzmann equa- 
tion, its exact form can be justified a posteriori within the 
framework of IDM (see appendix C). 

The general solution of equation (0 for the total energy- 
density, using eq. dTTb . is: 



p(a) = Cia'" + 



a F(a) 



[C2-J"x-if(x)F(x)dx] ' 
where the kernel function F(a) has the form: 



F(a) - exp 



-2kC 






H(x) 
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Note that kC\ has units of Gyr -1 , while m, C\ and C 2 are the 
corresponding constants of the problem. Obviously, eq.dT2l> can 
be seen as 



p(a) = p c (a) + p (a) 



(14) 



where p c - C\a m is the density corresponding to the residual 
"matter creation", resulting from a possible disequilibrium be- 
tween the particle creation and annihilation processes, while p 
can be viewed as the energy density of the self-interacting dark 
matter particles which are dominated by the annihilation pro- 
cesses. This can be easily understood if we set the constant C\ 
strictly equal to zero, implying that the creation term is negligi- 
ble, which reduces the current solution ("eqfT4b to that of eq.®. 
Note that near the present epoch as well as at late enough times 
(or » 1), as also in Case 2, the p evolves as the usual dark 
matter (see also Weinberg 2008). Finally, if both k and ¥ tend 
to zero, the above cosmological model reduces to the usual 
Einstein-deSitter model (Case 1). 

Note that, due to p > 0, the constant C 2 obeys the follow- 
ing restriction: 



C 2 > G(a) 



^ 



x^f(x)F(x)dx > 



(15) 
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Evaluating now eq.([T2t at the present time (a = 1, F(a) = 1), 
we obtain the present-time total cosmic density, which is: po = 
C\ + l/C 2 , with Ci > and C 2 > 0. 

Case 4: P - - y i'/3H: In this scenario we assume that the 
annihilation term is negligible [k - and f(a)=0] and the 
particle creation term dominates. Such a situation is mathe- 
matically equivalent to the gravitational DM particle creation 
process within the context of non-equilibrium thermodynam- 
ics Prigogine et al. (1989), the important cosmological reper- 
cussions of which have been studied in Lima et al. (2008 and 
references therein). Using our nomenclature and k — 0, eq.© 
becomes a first order linear differential equation, a general so- 
lution of which is: 



p(a) - a 



f 



x R(x)dx + Ci 



(16) 



The negative pressure can yield a late accelerated phase of the 
cosmic expansion (as in Lima et al. 2008), without the need 
of the required, in the "classical" cosmological models, dark 
energy. 

Below, we investigate the conditions under which eqs. (fT2l) 
and ( TTol l could provide accelerating solutions, similar to the 
usual dark energy case. 

4. Case 3:P = (kp 2 - W)/3H 

4.1. Conditions to have an inflection point and galaxy 
formation 

In order to have an inflection point at a - aj we must have 
&i — (see eq.[6]l. The latter equality implies thatp + 3P = 
should contain a real root in the interval: a e (0, 1). Therefore, 
with the aid of eq. lfTZl) . (0 and (fTTl) . we derive the following 
condition: 



a^FiH + 2kC\ a m ) kq^F 



-6172 



C 2 -G 



(C2-G) 2 



(m + 2)Cia"'H = 0, (17) 



from which we obtain that m > —2 (where C\ > 0, k > 
and C2 - G > 0). Evidently, if we parametrize the constant 
m according to m - -3(1 + widmX we obtain the condi- 
tion: widm < -1/3, which means that the current cosmolog- 
ical model can be viewed as a viable quintessence dark-energy 
look-alike, as far as the global dynamics is concerned. Indeed, 
we remind the reader that the same restriction holds for the 
usual dark energy model in which Pq - wpQ (w = const). 
Since the avenue by which the IDM model provides cosmic ac- 
celeration may appear slightly involved, we present in appendix 
D its correspondence to the usual dark energy models. 

Furthermore, in order to have growth of spatial density fluc- 
tuations, the effective DM part should be capable of cluster- 
ing and providing the formation of galaxies, while the effective 
dark energy term should be close to homogeneous. Indeed in 
our case the effective term that acts as dark energy is homoge- 
neous in the same sense as in the classical quintessence, while 
the Kp 1 term slightly modifies the pure DM evolution. In any 
case the interacting DM term after the inflection point tends to 
an evolution oc a 3 . During the galaxy formation epoch at high- 
z's we expect (due to the functional form of the DM term) that 
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Fig. 1. Left Panel: The Q2,o - kC\ solution space provided by 
fitting our model to the early-type galaxy Hubble relation of 
Simon et al (2005). Right Panel: The corresponding ^2,0 - M x 
solution space. 

the slope of the interacting DM term is not far from the clas- 
sical DM evolution (we will explore further these issues in a 
forthcoming paper). 

4.2. Relation to the Standard A Cosmology 

As an example, we will show that for m = (or widm = — 1) the 
global dynamics, provided by eq. (fT2l . is equivalent to that of 
the traditional A cosmology. To this end we use dt = da/{aH) 
and the basic kernel (eq.[T3l) becomes: 



F(a) - exp 



-2kC 



f 



1 



xH(x) 



-dx 



-2KC,(t-to) 



(18) 



where to is the present age of the universe. In addition, the in- 
tegral in equation ( fT2l see also eqIT5l> takes now the following 
form: G(a) = -kZU) andZ(f) = fVV 2 * Cl <'-">'. Note that at 
the present time we have G(l) = 0. Therefore, using the above 
formula, the global density evolution (eq. fT2l can be written: 



p(a) - C\ + a 



e -2KC X (t-t ) 

[Ci - G(a)] 



(19) 



As expected, at early enough times (t — > 0) the overall den- 
sity scales according to: p(a) oc cr 3 , while close to the present 
epoch the density evolves according to: 



P(a)-Ci + — , 
C2 



(20) 



which is approximately the corresponding evolution in the A 
cosmology, in which the term C\ acts as the constant-vacuum 
term (pa) and the I/C2 term acts like matter (p m ). 

Note that the effective pressure term (eq. [3J, for k — » 0, be- 
comes: ¥ ~ 3CiH, which implies that: P ¥/3i7 = -C\. 

Therefore, this case relates to the traditional A cosmology, 
since C\ corresponds to p\ (see eq. l20l . We now investigate 
in detail the dynamics of the m - model. 

From eq.dT9t. using the usual unit-less Q-like parameteri- 
zation, we have after some algebra that: 






Qi Q2,oa e 



3--2*CiC*-«d) 



Ql,0 + KCl£l2,oZ(f) 



(21) 



Spyros Basilakos 1 and Manolis Plionis 2 : Dark matter interactions and dark energy 



with Q lj0 = 8ttGCi/3HI and Q 2;0 = 87rG/3#JC 2 , which in 
the usual A cosmology they correspond to O a and Q. m , respec- 
tively. 

We can now attempt to compare the Hubble function of 
eq.(t2Tb to that corresponding to the usual A model. To this end, 
we use a x 1 minimization between the different models (our 
IDM eqflT] or the traditional ACDM model) and the Hubble 
relation derived directly from early type galaxies at high red- 
shifts (Simon, Verde, & Jimenez 2005). For the case of our IDM 
model we simultaneously fit the two free parameters of the 
model, ie., Q2,o and kC\ for a flat background (Oio = 1 — £2 2 ,o) 
with //o = 72 km/sec/Mpc and to = H^ ' ^ 13.6 Gyrs (roughly 
the age of the universe of the corresponding A cosmology). 
This procedure yields, as the best fitted parameters, the follow- 
ing: Q.2$ - 0.3*oog and log(/fCi) - -9.3 (with stringent up- 
per limit =b -3, but unconstrained towards lower values) with 
X 2 l&.f- - 1.29 (see left panel of Fig. 1). Using eq.(4) we can 
now relate the range of values of kC\ with the mass of the DM 
particle from which we obtain that: 



M x = 



1.205 x 10~ 1Z (cru) 



kC x 



10 



-22 



GeV, 



(22) 



(see also right panel of Fig. 1) and since kC\ is unbound to- 
wards small values, it is consistent with currently accepted 
lower bounds of M x (~ lOGeV) (eg., Cirelli et al. 2009 and 
references therein). The corresponding Hubble relation (Fig. 
2), provided by the best fitted free parameters, is indistinguish- 
able from that of the traditional ACDM model, due to the very 
small value of kC\ - 10~ 9 3 . For completion we also show, 
as the dashed line, the IDM solution provided by M x ~ leV 
{kC\ — 10~ 3 ), which is the stringent lower bound found by our 
analysis. In this case the predicted Hubble expansion deviates 
significantly from the traditional A model at small a values in- 
dicating that it would probably create significant alterations of 
the standard BBN (eg. Iocco et al. 2009 and references therein). 
Although the present analysis does not provide any impor- 
tant constraints on M x (within our model), we plan to use in 
the future a large number of cosmologically relevant data to at- 
tempt to place stronger M x constraints, also for the general case 
of section 4.1. 



5. Case 4: P = -Y/3# 

In this section we will prove that for k — (negligible anni- 
hilation), the global dynamics resembles that of the traditional 
quintessence cosmology (w =constant). Indeed, using again the 
phenomenologically selected form of *P, provided by eq. dTTb . 
we have R(a) - C\(m + 3)a m ~ 1 , It is then straightforward to 
obtain the density evolution from eq.([T6li. as: 



p(a) 



Da- 3 +Cia" 



(23) 



where D - Ci - C\ . The conditions under which the current 
model acts as a quintessence cosmology, are: T> > 0, C\ > 
and widm = — 1 - m/3, which implies that in order to have an 
inflection point, the following should be satisfied: wjdm < -1/3 
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Fig. 2. Comparison of the Hubble function provided by the tra- 
ditional ACDM model, which coincides with our ra = model 
(for the best fit of the two free parameters - see text). The 
dashed line corresponds to our ra = IDM model for the high- 
est kC\ bound, provided by our fitting procedure (~ 10~ 3 ). The 
dot-dashed line corresponds to our k — IDM model (Case 
4) for the best fitted parameters (m =* -0.3 and Q2,o - 0.28). 
Finally, the points correspond to the observational data (Simon 
et al. 2005). 



or m > -2 (see appendix D). Notice, that the Hubble flow is 
now given by: 



H . 

— = Q 2 ,oa 

tl{) 



+ fil,0tt" 



(24) 



with O 2 ,o = 87TGD/3H 2 and Q 1)0 = %nGCil3H*. Finally, 
by minimizing the corresponding x 2 ( a s in section 4.2), we 
find that the best fit values are Q 2 ,o - 0.28 and m - -0.30 
(widm — -0.90) with ,y 2 /d.f. = 1.29. The corresponding 
Hubble flow curve is shown in Fig. 2 as the dot-dashed line. 
Note that this solution is mathematically equivalent to that of 
the gravitational matter creation model of Lima et al. (2008). 

6. Conclusions 

In this work we investigate the evolution of the global density 
of the universe in the framework of an interacting DM scenario 
by solving analytically the collisional Boltzmann equation in 
an expanding Universe. A disequilibrium between the DM par- 
ticle creation and annihilation processes, whichever its cause 
and in favor of the particle creation term, can create an effective 
source term with negative pressure which, acting like dark en- 
ergy, provides an accelerated expansion phase of the universe. 
There are also solutions for which the present time is located af- 
ter the inflection point. Finally, comparing the observed Hubble 
function of a few high-redshift elliptical galaxies with that pre- 
dicted by our simplest IDM model (m = 0), we find that the 
effective annihilation term is quite small. In a forthcoming pa- 
per we will use a multitude of cosmologically relevant obser- 
vations to jointly fit the predicted, by our generic IDM model, 
Hubble relation and thus possibly provide more stringent con- 
straints to the free parameters of the models. We plan also to 
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derive the perturbation growth factor in order to study structure 
formation within the IDM model. 
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Appendix A:The effect of the decay products 

Here we attempt to investigate in the matter-dominated era, whether 
the possible radiation products due to dark matter interactions can af- 
fect the global dynamics. A general coupling can be viewed by the 
continuity equations of interacting dark matter Pidm and residual radi- 
ation 6p r , 



Dark matter interactions and dark energy 

- Case 2: For the case where: 

R(x) = h'(x) with g(x) = -f(x)h(x) 



(32) 



the particular solution is h(x) [in our case we have h(x) 
-3k~ ] H(x)]. The general solution now becomes: 



y(x) = h(x) + O(jc) 

where 

O(.y) = exp 



C 2 



- I f(u)<t>(u)du 



I f{u)h(u)du 



(33) 



(34) 



+ 3H(t)p lDM + Kp( DM - ¥ = Q 



At 



dSp r 

dt 



+ 4H(t)6p r = -Q 



(25) 



(26) 



where Q is the rate of energy density transfer. If Q < then the 
IDM fluid transfers to residual radiation. As an example we can use 
a generic model with Q = -eSp r , where e > 0. Thus, eq.J26b has an 
exact solution 



In this framework, using eq. (8) we finally get ^(x) = 
xH(x)R(x) = -3k- [ xH(x)H'(x). 

Note that the solution of Case 1 is the only one providing a oc a~ 3 
dependence of the scale factor (see eqs. ll2l|19| andl20l. 



Appendix C: Justification of functional form of ¥ 

Suppose that we have a non-perfect cosmic fluid in a disequilibrium 
phase with energy density p. Then from the collisional Boltzmann 
equation, we have: 



Sp r = Spaa e 



■4 -,*-'<>) 



(27) 



¥ = p + 3Hp + Kp 1 = —aH + 3Hp + Kp 2 
da 



(35) 



with t the present age of the Universe. This shows that the contribu- 
tion of the residual radiation to the global dynamics is negligible in the 
past, since there is not only the usual oc or 4 dependence of the back- 
ground radiation but also a further exponential drop, and thus 2-0 
Therefore we conclude that we can approximate the total energy- 
density with that of the interacting dark-matter density (p = Pidm)- 
Note, that 1/e can be viewed as the mean life time of the residual 
radiation particles. 

Appendix B: Solutions of the Riccati equation 

With the aid of the differential equation theory we present solutions 
that are relevant to our eq.([7j- In general a Riccati differential equation 
is given by 



Furthermore, we assume that for a convenient period of time the cos- 
mic fluid, in an expanding Universe, is slowly diluted according to 
p ~ C\a m (m < 0). From a mathematical point of view, the latter as- 
sumption simply means that a solution of the form oc a" 1 is a particular 
solution of the Boltzmann equation. Therefore, we have finally that: 



V^C\(m + 3)a m H + KC]a 2 "' 



(36) 



Appendix D: Correspondence between our model 
and the usual Dark Energy models 

We remind the reader that for homogeneous and isotropic flat cos- 
mologies (fi ra + Q.q = 1), driven by non relativistic DM and a DE 
with a constant equation of state parameter (w), the density evolution 
of this cosmic fluid can be written as: 



y = f(x)y 2 + g(x)y + R(x) 



(28) 



p(a) = p^oa-' + p Q ,oa- 3{i+w) 



(37) 



and it is fully solvable only when a particular solution is known. Below 
we present two cases in which analytical solutions are possible. 



- Case 1: For the case where: 



R(x) = Cm*"- 1 - C 2 l x 2m fM ~ C^'gix) 



(29) 



the particular solution is x"' and thus the corresponding general 
solution can be written as: 



I f{u)<b{u)du 



y(.x) = CiX m +®{x) 

where 

*(x) = exp J (2Ci u'"f(u) + g(u)) du 



(30) 



(31) 



and C{,Ci are the integration constants. Using now eq.(8) we get 
T(x) = xH(x)R(x) = C\{m + 3)x"'H{x) + kCW" 1 . 



where p m0 and p G ,o are the present-day DM and DE densities, respec- 
tively. 

The necessary criteria to have cosmic acceleration and an inflec- 
tion point in our past (/, < t ), are: (a) P < and (b) a = 0, which 
leads to the conditions: 

- Dark Energy models: P = P m + Pq = wp Q < 0, P m = with 

w<-l/3. 

- IDM models: P = Kp 2 - "P < and m > -2 (or widm < -1/3 see 
section 4). 
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